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INTERACTIONS  BETWEEN  OCEAN 
SURFACE  WAVES  AND  CURRENTS 


1.  INTRODUCTION 

A  surface  wave  field  may  be  significantly  altered  by  interaction  with  currems.  These  currents  can 
be  induced  by  many  mechanisms  in  the  ocean.  A  potentially  important  wave-current  interaction  results 
from  the  interaction  of  ocean  surface  gravity  waves  with  the  wake  produced  by  a  moving  object  or 
vessel.  Visual  observations  suggest  that  in  the  wake  belund  a  surface  ship  there  exists  a  'dead-water' 
region  which  is  characterized  by  the  relative  absence  of  short  wavelength  waves.  Tlus  r^on  often 
extends  to  very  large  distances  teh’nd  the  body.  Relative  to  the  ambient  sea,  the  wake  region  appears 
to  be  very  smooth.  This  is  further  displayed  by  some  of  the  synthetic  aperture  radar  (SAR)  images 
made  of  surface  ship  wakes. 

In  this  report  we  will  evaluate  a  potential  model  to  predict  the  far-field  development  of  this  quies¬ 
cent  r^on.  We  will  consider  the  possibility  that  it  is  the  interaction  of  the  ambient  and  Kelvin  wave 
fields  with  the  momentum  wake  of  the  vessel  which  is  responsible  for  the  generation  of  this  region.  In 
cct^ction  with  the  conceptual  model  of  the  total  wake  flow  field  previously  developed  by  Skop  [1], 
and  by  Cooper  and  Skop  [2]  we  will  assume  that  the  leading  order  approximation  to  this  interaction 
problem  is  obtained  from  the  interaction  of  the  turbuletit  momentum  wake  as  superimposed  on  an  oth¬ 
erwise  quiescent  background  and  the  actual  background  wave  field.  The  results  of  Inis  phase  of  the 
study  will  provide  some  initial  estimates  of  the  effectiveness  of  this  mechanism  for  the  generation  of 
the  dead  water  region.  Since  our  interest  is  in  the  short  waves,  we  will  further  assume  that  the  wave- 
field  of  interest  is  composed  towlly  of  short  wavelength  waves.  The  deep  water  dispersion  relation  may 
be  utilized  throughout  and  the  vertical  variatioirs  in  wake  velodty  may  be  neglected.  The  surface  veloc¬ 
ity  field  produced  by  the  wake  may  then  be  considered  as  a  current  field  in  its  interaction  with  the  sur¬ 
face  waves. 

Some  of  the  earliest  work  on  the  problem  of  the  interaction  of  regular,  linear,  surface  waves  and 
currents  is  that  of  Longuet-Higgins  and  Stewart  [3,41.  This  work  as  well  as  other  work  on  regular  wave 
interactions  is  summarized  in  the  comprehensive  review  article  by  Peregrine  [51.  Peregrine  [6]  also 
developed  a  model  to  estimate  the  overall  effect  of  the  wake  velocity  field  of  a  ship  upon  the  Kelvin 
waves  generated  at  the  ship’s  stem  and  was  able  to  demonstrate  that  for  a  sufficiently  strong  'effective* 
wake,  the  envelope  of  the  stem  wave  is  signific^fly  reduced  from  its  ideal  no:,  interactive  value.  More 
recently.  Peregrine  and  Thomas  {7]  have  considered  the  interaction  for  finue-ampUtude  symmetric 
waves.  The  finite-amplitude  wave  properties  are  computed  using  the  method  of  Longuet-Higgins  18] 
for  predicting  the  energy  and  wave  action  densities  and  fluxes.  Studies  of  irregular  or  random  wave 
field  interactions  were  conducted  by  rrumerous  workers  (9,  10,  11,  12,  13].  The  work  of  Tayfun  et  al. 
ilO]  and  James  el  al.  ill]  are  most  relevant  here  as  they  form  the  basis  fur  the  work  to  be  presented  in 
this  report 


Maauscnpl  tptroved  Januaiy  20,  1986. 


II.  ANALYSIS 
A.  Wake  Velodty  Field 

Far  downstream  of  the  trailing  edge  of  a  moving  object  we  expect  the  wake  to  possess  only  hor¬ 
izontal  velocity  components  u  and  v.  The  vertical  velocity  component  tv  -  0  and  the  variations  of  a 
and  V  will  be  predominantly  in  the  transverse  direction,  y.  Very  slow  variation  in  the  downsueam  coor¬ 
dinate  X  is  expected  relative  to  the  y  variation.  We  therefore  consider  the  fundamental  problem  a  wake 
velocity  field  given  by 

V-  (u(y).  vfvl  0) 

where  y  is  transverse  to  the  direction  of  motion  of  the  body,  —x.  If  V  is  to  be  a  legitimate 
incompressible  wake  velocity  field,  then  V  •  V  -  0  from  continuity,  implies  that  v(y)  =  0*  or 
V  -  (u(y),  0,  0).  That  is,  far  downstream  of  the  trailing  edge  the  wake  is  expected  to  be  described  by 
parabolic  equations  with  a  single  velocity  component  u(y)  being  predominanL  However,  if  the  body  is 
not  infinite  in  the  vertical  direction  or,  equivalently,  if  there  is  a  finite  velocity  component  ia  the 
vertical  direction  then  a  legitimate  wake  field  is  V  -  (i/(y),  v(y),  >v(y,  z)).  Then,  firm  csuimuity, 

__  .im 

By  Bz' 


It  is  expected  funner,  Garrett  ti4i,  that  in  the  lar  wake,  v/u  <<i.  Therefore,  the  simplest 
interaction  wake  field  and  one  which  appears  to  bear  close  connection  to  observation  is  a  current  profile 
V  -  u(y).  The  effects  of  a  finite  transverse  velocity  component,  v,  will  also  be  addressed  in  some  of 
this  work. 

B.  Ware  Kinematic  Reiations 

In  this  section  we  will  derive  the  kinematic  relations  which  are  satisfied  by  the  individual  Fourier 
modes  of  the  wave  field  spectrum  in  their  interaction  with  the  surface  current.  Each  mode  is  approxi¬ 
mated  as  a  nearly  plane,  slowly  varying,  small  amplitude  wave  field.  Subsequently,  the  results  of  this 
section  will  be  utilized  to  develop  the  relevant  equations  to  describe  continuous  tandom  wave  fielos 
which  are  appropriate  for  real  ocean  waves 

Consider  a  nearly  plane,  slowly  varying  wave  field  This  field  may  be  approximately  described  by 
an  amplitude  function  A  (x,  r)  and  a  phase  function  Six,  r).  These  two  functions  satisfy 

ill  -  A{x,  r)  exp  /(Six,  r)) . 

In  analogy  with  plane  waves  where  S  s  (k  -  x  -  wr),  we  draw  the  following  definitions  for  wave 
number  and  frequency 

k  -  VS  (1) 

and  <fl  —  —  .  (2) 

From  Eq.  (1)  we  have  V  x  k  -  0  and  by  eliminating  S  from  (I)  and  (2)  we  obtain 

^  +  Voi-0.  (3) 

Bt 


sate  of  v  ts  d«termin«d  from  the  slow  variaiion  9u/dx  end  cs  neglisibte  to  lowest  order 


TUs  is  the  well  known  conservation  of  crests  equation.  To  complete  the  description,  a  disperaon  rela¬ 
tion  is  required  which  also  relates  a  and  k.  For  traveling  waves  which  propagate  relative  to  a  coordi¬ 
nate  system  which  is  fixed  in  the  fluid  we  have 

a  <rQc).  (4) 

Here  a-  denotes  the  intrins’c  frequency.  If  the  medium  possesses  a  slowly  varying  current  U(x,  t)  then 
with  respect  to  a  stationary  coordinate  system  we  have 

<a-<r(k)  +  k-u(i,r).  (5) 


If  we  take  the  time  derivative  of  Eq.  (S)  we  obtain 

or  from  Eq.  (3)  we  have 


dr  3k  3t  or  dr 


"«>  .  I  do- 


dr  3k 


■,T+  ^  +  «  Vw-k-^ 


3U 


dr 


do- 


where  is  recognized  as  the  group  velocity.  That  is, 
dk 


£q-cation  (7)  is  expressed  as 


c 

*  3k  • 


dr  *  dr 


(6) 


(7) 


(8) 

(9) 


or 


'I 


(10) 


where  represents  a  time  derivative  along  rays.  Equation  (10)  therefore  indicates  that  the  total  fre- 
at 

quency  a  is  conserved  along  rays  for  a  steady  current  or  when  observed  in  a  coordinate  system  where 
the  current  is  steady.  The  wake  field  is  steady  when  referred  to  a  body-flxed  coordinate  system.  The 
transformation  to  a  body-flxed  coordinate  system  is  therefore  advantageous  since  the  total  frequency  a 
of  the  waves  will  then  be  invariant  during  the  interaction. 


The  ambient  wave  field  also  must  be  characterized  with  respect  to  the  same  body-flxed  coordinate 
system.  Consider  an  ambient  traveling  wave  at  a  frequency  u.  and  wave  number  k.  where 

«>»  “  <T(k»)  (11) 

If  the  body  which  is  generating  the  wake  is  translating  at  a  uniform  velocity  then  vrith  respect  to  a 
body-flxed  observer  the  frequency  would  be  given  as 

01 "  oiw  —  k,»  *  Do 

«  -  o-(*„)  -  k.  ■  (12) 


and  the  waxe  velocity  field  U  —  — Uj  +  d  where  d  is  the  wake  field  relative  to  a  fixed  observer.  In 
Equation  (12)  we  have  used  the  fact  that  k.  is  unchanged  in  the  Galilean  transformation.  The  fre¬ 
quency  0,  given  by  (12)  is  the  total  frequency  which  is  invariant  through  the  interaction.  That  is,  by 
combining  (12)  and  (S)  for  a  steady  wake  field  we  obtain 

o:  -  cr(k.)  -  k»  •  li»  -  o(k)  -b  k  ■  U.  (13) 


3 


(14) 


If  we  sp^alize  to  lin^  dwp  water  waves  then  we  may.write 

<r(k) -(*?)« 

and  Eq.  (13)  becomes 

(*:»«)’'?-*„, £4 -  (fe)«  +  k  U- (*y)«-ilil/j+k-fi  (15) 

where  we  have  taken  in  the  negative  x  direction. 


k. 


By  combining  Eqs.  (3)  and  (5)  and  eliminating  a  we  obtain  the  ray  equation  for  the  wave  number 


(16) 


It  is  informative  to  examine  the  two  component  equations  of  Eq.  (16).  In  the  x  direction  we  have 


—  +  (u  +  t,)  —  +  (v-*-  c;.)  —  -  -k.  iil  -  t.  ^ 
a,  ^  ^  dx  3y  ^  9x  ^  a* 


9a 


and  in  the  y  direction 


9x 


By 


93- 


9y' 


(16a) 


(161,) 


From  Eq.  (16a)  we  note  that  for  £/-  (u(y),  v(y),  w),  the  right  ade  va.oishes  and  we  have 

dk, 

-  0  or  is  conserved.  The  variation  of  is  described  by  Eq.  (16b).  Tht  f  ardmate  y  is  a  hor- 
iiontal  direction  transverse  to  the  wake  axis  (x). 

For  kj  -  constant,  the  terms  involving  l/j  cancel  from  both  sides  of  Eq.  (IS)  and  wc  find  that 
Eq.  (IS)  provides 

(k»s)'^^- (fcs)’'*  +  k-5.  (17) 


This  equation  indicates  that  the  vave  refraction  is  governed  by  the  viscous  or  turbulent  velocity  for  the 
case  where  Ay  —  constant  Tha  idnienatics  of  the  refraction  process  are  therefore  governed  by 

(kcg)'*’^  “  +  Ay  C  +  ^v  (18) 

and 


Ay  "•  ■-  constant 


(19) 


From  Eq.  (IS)  ae  have 


Equation  (IS)  provides 

(A.g)'/^- 


Ay  *  “  /;.» 

Ay  ••  A  cos  9  — 


err,  -•  A'sin  9 
sin  9,,  cos  9 
sin  9 


Ar_  sin  9. 


sin  9 


+  5  A^  sin  9„  —  V 


A,  sin  9» 
sin9 


cos9 


(20) 


(21) 
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where  0  represents  the  angle  between  k  and  the  y-axis  and  for  v  >  0  we  consider  only  <  0  as  indi¬ 
cated  by  the  minus  sign  preceding  the  v  terms.  Therefore,  — y  <  0  <  y,  *,  <  0  and  Eq.  (21)  may 
be  written  as 

IW  I  .1 


+  sin0_ia* 


.  y.cos0 
sin0 


If  we  define 


f  =  sin0 
*  ~  sin  0« 


then  Eq.  (22)  may  be  expressed  as 

-  sin  0„u*)  -  {  -  V*  -»/l-f*sin0».  (25) 

We  note  that  in  the  far  field  where  u*  and  v*  vanish,  sin  0  — •  sin  0«  or  {  —  1.  Also,  in  the  limit 
where  V*  —  0  we  recover  Eq.  (2.151)  of  Peregrine  151  which  is  consistent  with  v  —  0. 

Equation  (25)  is  valid  for  refraction  of  the  wave  field  by  a  general  two  component  current 
(u(y),  v(y)).  In  what  follows,  we  will  neglect  v  in  comparison  to  u.  In  future  work  on  tlus  problem 
the  effects  of  finite  v  variation  may  be  addressed.  The  basic  problem  described  by  Eq.  (25)  with  v  —  0 
is  the  interaction  of  a  wave  field  with  a  parallel  surface  current  shear  u  (y). 


The  variation  of  wave  amplitude  due  to  the  refraction  may  be  determined  by  the  conservation  of 
wave  action.  For  the  steady  case  considered  here  this  is  expressed  as 


V|(5-!-cP'^J-0,  (26) 

where  ^  represents  the  wave  energy  density  for  linear  waves  of  amplitude  a,  and  o-  once 

again  represents  the  intrinsic  frequency  relative  to  the  fluid.  For  the  parallel  shear  flow  considered, 
u(y),  Eq.  (26)  implies  that 


Cj,  yj  »  constant 


cos  0  yj  —  c,  cos  0  y  ^^--j  -  ( 


^  1 

For  linear  deep  water  waves  —  -  -r-,  so  that 
O’  2k 


and  therefore 


CpF  _  COS0  pga^  sin 0  cos  0pgo^ 
O’  4  *r  4*, 


a^pg  sin  20  . . 

^  P  *  vOnst  ttit« 


5 


It  am  te  t^t  the  wave  amplitude  becomes  unbounded  when  sin  29  0  or  0  ±-~,  0.  As  indi- 

cated  by  Peregrine  [S],  0  —  ±  -y  correspond  to  a  caustic  where  the  plane-wave  assumption  breala 
down.  For  0-0.  wave  steepness  becomes  unbounded. 

C.  Random  Wave  FieU 

The  above  results  and  discusaon  have  been  associated  with  the  behavior  cf  the  indivddual  Fourier 
components  of  a  wave  field.  In  practical  problems,  however,  we  are  usually  confronted  by  a  continuous 
spectrum  of  random  waves.  The  refraction  of  wave  fields  of  this  type  has  been  the  sulgect  of  a  great 
deal  of  recent  work  (9,  10,  11,  12,  131.  These  random  fields  are  typically  described  by  the  inhomo¬ 
geneous  wave  number  spectral  density  d>(k;  x)  which  is  defined  by 

dElpg  s  d'ik:  x)  d’k  -  i(t(k;  x)  *,<*2  -  ^kdkda  (28) 

where  dE  represents  the  wave  energy  at  wave  number  k.  The  total  energy  associated  with  the  local 
mean  square  value  of  the  free  surface  displacement  is  obtained  by  int^rating  dE  over  ail  of  k  space. 
That  is. 


£  -  /  /  di(k,  x)  dkfdki  -  <{’(x)> . 


where  £  represents  the  mean  square  surface  elevation  and  is  defined  by  £  —  £/pg.  But  the  energy  can 
also  be  defined  in  terms  of  a  frequency-direction  spectral  density  <b(ti>,  0,  x)  as 


£  —  J*  J"  4>(u,  0;  x)  &>  du  dB. 


Therefore  from  Eqs.  (29)  and  (30) 

*(k,x)-^-^q.(w,0;x).  (31) 

k  dk 

As  previously  demonstrated  for  steady  wave-current  interactions,  a  is  constant  and  therefore  the  fre¬ 
quency  direction  spectral  density  <l>(<a,  0;  x)  is  the  more  convenient  variable. 

From  Eq.  (28)  we  can  derive  the  consequences  of  conservation  of  wave  action  (for  each  wavelet). 
That  is  for 

dE/pg  -  ili(k,  x)  d’k 

the  wave  action  is  written  as  (’ftd^Vtr  and  therefore  conservation  of  wave  action  demands  that 


^  ^d^  +-^  (U,  +  c^)  -4-  a]  -  0 

Or  <r  9xi  a- 


It  can  be  shown,  Le  Blond  and  Mysak  IS]  and  James  et  al.  Ill],  that 

(</=k)  +  (((/,  +  <i,)  d^)  -  0.  (33) 

at  dx,  ^ 

Therefore,  combining  this  result  with  Eq.  (32)  we  obtain  a  conservation  law  for  the  spectral  action  den¬ 
sity  which  is  given  as 


|:N  +  (<i,-n/.)^kl-o 


Equation  (34)  indicates  that  the  spectral  density  of  wave  action  ^/<r  is  conser.'ed-along>tays.  The 
reader  should  note  the  difference  between  Eqs.  (32)  and  (34).  Por  steady  interactions  where  a  is  also 
conserved  along  rays  it  is  much  more  convenient  to  study  die  spectral  modifications  of  •!>,  the 
frequency-direction  spectral  density.  From  Eq.  (31)  we  End  that 

<r  fe  <ft  <r  'W 

Therefore  for  steady  interactions  where  a  and  are  conserved  along  rays,  Eq.  (3S)  implies  that 

O  ■' 
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r — - IS  conserved. 

k  dk  <r 

These  properUes  were  utilized  by  James  et  al.  (Ill  in  their  studies  of  the  refraction  of  a  random 
wave  field  by  a  parallel  surface  shear  current  They  included  the  effects  of  wave  breaking  and  wave 
reflection  and  generated  the  resulting  mean  square  surface  displacement  given  by  Eqs.  (29)  or  (30). 
Thdr  computations  were  for  some  rather  drastic  and  hence  non-physical  assumptions  relating  to  the 
ambient  frequency-direction  spectral  density.  These  assumptions  simplified  the  computations.  It  would 
appear  however,  that  it  is  possible  to  utilize  more  realistic  ambient  spectral  densities  within  a  similarly- 
ba^  formulation  of  the  problem.  Furthermore,  it  is  possible  to  extract  specific  information  about  tte 
surface  energy  contributions  from  particular  wavelength  ranges  in  the  interaction  me  and  thereby 
be^n  to  address  the  question  of  the  origins  of  the  dead  water  re^on.  Most  importantly,  't  is  posable  to 
evaluate  from  the  spectral  density  modifications,  free  surface  spatial  realizations  cortespo  ‘ding  to  vari¬ 
ous  random  sea  state  models. 

The  techniques  developed  by  James  et  al.  (11]  enable  the  calculation  of  the  modification  of  the 
energy  associated  with  the  mean-square  free  surface  elevation  of  the  incoming  wave  field  due  to  the 
refraction  caused  by  a  surface  current.  They  assumed  a  simplified  incoming  spectral  density  distribu¬ 
tion  which  was  taken  to  be  independent  of  frequency  and  direction.  That  is,  they  assumed 

t-fw..®-)  -  1. 

This  approximates  the  incoming  wave  spectrum  as  one  similar  to  that  of  whits  noise.  In  addition,  all 
the  kinematics  of  the  refraction  process  were  referred  to  the  initial  plane  (u.,P_)  where  the  incoming 
energy  for  a  given  wavelet  is  proportionai  to  the  polar  area  in  that  ph  ' 

dEm  ~  ~  (36) 

when  <I>„  -  1.  From  Eqs.  (34)  and  (35)  we  can  express  the  ratio  of  the  local  frequency-direction  spec¬ 
tral  density  <I>  to  its  original  undisturbed  value  <1>.  along  a  ray.  That  ratio  is  given  by 


.  _  <P(m,9;x) 

“  <I>»(w-,«.;x>) 


where  x  and  x_  lie  on  the  same  wave  ray.  From  Eq.  (18),  we  have 

cii  —  <rik}  +  *,«;  —  cr„.  ( 

For  linear  deep  water  waves,  <r(.k)  -  (kg)*^^ .  Therefore,  combining  Eqs.  (38)  and  (37)  we  obtain 
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or,  since  -.0  by  de{milio%  vre  h»re 


A 


Furthermore,  for  the  case  where  the  transverse  component  of  2,  is'negligible  (2  ••  0)  we  find  that 

ts 


I  tt)  [  Cl) 


(41) 


Following  lames  et  al.  [11]  we  adopt  the  following  normalization  for  k  and  oi  based  upon  the  maximum 
current, 


and  we  define  a  velocity  scale  A  by 


A  =  A(y)  - 


As  a  consequence,  we  find  from  Eqs.  (20)  and  (25)  that 


*i“»-.^.a_Ae..sinO«. 

t  * 


sin  d  Am 

Defining  ij  —  u  *  sin  0  and  t).  ~  <u  *  sm  0.  we  have 


(42) 


JU.. 


r;-  (1  —  Asjm)* 


(43) 


which  expresses  the  local  wave  direction  in  terms  of  the  initial  direction.  We  can  also  obtain  tlie 
inverse  solution  by  solving  the  quadratic  for  t)..  This  yields 

where  we  have  selected  that  root  which  provides  i]“’  iJm  at  A  —  0  . 


1  +  2A')}  —  -Jl  +  4Ai7 


(44) 


In  terms  of  the  above  normalization,  Eq.  (41)  for  the  frequency-direction  spectral  denaty  ratio 
becomes 


<p  ,  (1  —  AiJm)^ 

C"  "  (1  +  A„_)  ’ 

and  if  we  substitute  the  inverse  refraction  Eq.  (44)  this  ratio  then  becomes 

.  _  [vi-h4ATi-  ll*  1 
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(45) 

(46) 


Equations  (45)  and  (46)  are  extremely  useful  for  defining  the  Fourier  decomposition  of  the  surface 
elevation  with  respect  to  the  initial  end  the  physical  frequency-direction  planes  respectively. 


D.  Fourier  Spectrum 

As  discussed  above,  a  random  wave  field  may  be  characterized  by  means  of  the  spectral  density. 
To  describe  the  refraction  of  these  waves  by  a  steady  current,  the  frequency-direction  spectral  den«ty  is 
most  convenient  since  the  total  frequency  is  conserved  in  such  an  interaction  as  demonstrated  above. 
From  Eq.  (30)  we  have 

dE  —  <t>((i>,d;z)udui/9 

where  dE  represents  the  mean  square  surface  displacement  energy  contribution  for  the  region  between 
9  and  9  +  do  and  <a  and  a  -b  da.  The  tohil  energy  is  obtained  by  integration  over  all  directions 
-*r  <  6  <  JT  and  over  all  frequencies  0  <  u  <  «.  That  is,  the  mean-square  free  surface  elevation  is 

d9  f47a) 

James  et  al.  [11]  choose  to  define  the  interaction  in  the  initial  frequency-direction  space  (a>_,  9j>  and 
therefore  compute  <J^>  as 

-  X.  X.  ^  W7b) 

employing  Eq.  (4S)  for  A .  In  this  work  we  will,  at  times,  find  it  more  convenient  to  utilize  Eq.  (47a) 
for  <{^>  and  to  evaluate  A  by  Eq.  (46),  thereby  evaluating  the  energy  integral  in  the  physical  plane 
(w,  9).  Eq.  (47b)  will  be  utilized  to  calculate  of  the  mean-square  surface  elevation  energy  wMle  we  trill 
find  it  more  convenient  to  utilize  Eq.  (47a)  in  physical  space  to  evaluate  the  Fourier  decomposition  of 
the  surface  elevation  to  be  discussed  below. 


It  therefore  follows  that  the  contribution  to  surface  elevation  corresponding  to  frequencies  of  Au 
about  aiy,  and  directions  A9  about  9,  is  given  by  (9, 17, 18,  20] 


(48a) 


“  •\J2^(.aj,9i,r)  aj  Ao  A9  • 

cos  Or  sin  9(  +  y  cos  9,)  —ajt  + 

From  Eqs.  (43)  and  (44)  we  have  k  -  kj(.k^k)  —  ** 


VI  +  4An  —  1 


2Ar| 


(48b) 


where 

7)  “  <uy  sin  9(. 

and  <l>u  represents  a  random  phase  shift  uniformiy  distributed  between  0  and  2ir. 

The  above  expressions  indicate  that  the  Fourier  coefficients  for  the  surface  displacement  can  be 
determined  from  the  local  value  of  the  frequency-direction  spectral  density  <I>(u,9).  The  two  expres¬ 
sions  derived  for  A,  Eqs.  (45)  and  (46),  representing  the  ratio  of  spectial  density  along  a  ray  can  be 
utilized  to  construct  the  local  spectral  density  from  its  initial  value  due  to  the  interaction  with  a  current 
distnbution.  The  surface  energy  can  be  obtained  by  means  of  integration  in  either  initial  space  (45)  or 
physical  space  (46).  Hence,  from  Eq.  (46)  we  have  in  physical  space 
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(49) 
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«'^-K TT^l^l 

cos  [k,jKx  sm  9,  +  }  cos  0,)  -  oj,t  -*-  6,j] 
where  6„  (w*,  6)  is  dc  ermined  from  Equation  ‘44>  as 


and  the  fuJl  free  surface  realization  or  displacement  at  a  given  time 

-  A.  t)  iSV 

In  Eq  (51),  v,e  have  c'.aU.a.ed  ail  C..  at  /  =  0  Equaiions  i50i  and  <511  form  the  basis  of  a  numerical 
procedure  m  ph>sical  space  to  determine  the  Surface  elevation 

V'ariou«  criteria  can  be  utilized  for  ihe  inclusion  of  Aave  breaking  m  the  model  The  one  selected 
by  James  ei  al  [111  represents  a  maximum  condition  The>  define  the  wave  breaking  boundary  as  the 
condition  that  the  local  component  of  the  rurren'  in  the  direction  u'’  the  wave  motion  is  equa'  but 
opposite  to  the  local  wave  group  veloci»y  in  dimension'ess  terms  this  is  equivale»'*t  to 

=  (S2a) 

or,  in  the  pb>  steal  plane 


Clearly,  these  cor.du.ons  are  conservative  they  are  independent  of  the  precise  value  of  the  wave 
ampiitdde  Tr^e  wave  amplitude  defined  by  these  conditions  is  infinite  In  a  continuous  refraction  pro* 
cess,  waves  break  we'l  before  these  boundaries  are  reached  Therefore,  this  represents  a  conservative 
condition  in  terms  of  the  unbroken  wave  energy  More  energy  is  lost  to  wave  breaking  than  is  indi¬ 
cated  by  these  conditions  It  may  be  more  realistic  to  establish  the  wave  breaking  boundary  based  on 
the  local  wave  steepness  (191  This  represents  a  entena  which  depends  both  on  the  initial  spectral  dis¬ 
tribution  and  the  integrated  growth  of  the  particular  spectral  component.  For  simplicity  in  this  paper  v/e 
will  adopt  the  conservative  criterion  established  above  However,  wc  recognize  the  limitations  of  the 
condition  and  will  reconsider  this  at  a  future  time 

The  wave-breaking  boundary  is  utilized  to  establish  limits  of  integration  for  mean  square  surface 
elevation  energy  <C^>  or,  equivalently,  limits  for  the  contributions  to  the  Fourier  coefficients  for  the 
surface  elevation.  Utilizing  Equation  (52b)  as  the  definition  of  the  wave  breaking  boundary  (9*)  we' 
obtain 


This  establishes  the  lower  integration  limit  for  9.  For  waves  with  components  in  the  direction  of  the 
current,  9  >  0,  interaction  with  the  current  can  produce  reflected  waves  Waves  which  are  about  to 
reflect  are  characterized  in  the  physical  plane  at  0  »  w/2.  In  the  initial  plane,  these  were  characterized 
by  lames  ct  al.  as  —  6*,,.  Therefore  the  9  iimiis  for  incoming  waves  in  Eqs  (49),  (50)  and  (51)  are 
given  as  0,^,^  <  ^  <  ^/2  0nan  the  minimum  of  9^  and  sin"*(-l/(l  +  Aw*)^)  obtained  from  equa¬ 


tion  (42)  al  9oc 


‘  The  outgoing  waves  can  be  shown  to  ooey  the  same  limits.  These  are  waves 


originating  on  the  other  side  of  the  wake  and  those  which  have  been  reflected  from  regions  beyond  the 


'  \r 


poini  in  quesuon  These  results  ma>  be  utilized  lu  generate  free  surface  realizations  within  the  interac- 
:.:"i  zone  ru'thcrmcrc,  it  is  possible  to  evaluate  other  surface  nroperues  which  may  be  important  in 
inierpreiirg  observations  For  example^  the  surface  slope  distribution  which  is  important  in  the  remote 
sensing  of  ocean  sur^’aces  is  expressible  Irom  a  knowledge  of  the  surface  elevation  Fourier  decompose 


Included  here  has  been  an  outline  of  procedures  to  obtain  the  spectral  density  distribution  within 
t’te  interaction  zone  of  a  wave  field  and  a  current  distribution  and  thereby  to  obtain  the  Fourier  spec¬ 
trum  and  a  spaMa>  realization  of  the  surface  elevation  This  work  is  currently  underway  and  further 
results  'vill  be  c'  ailable  shortly  both  for  the  simplified,  <I>»  *  and  for  more  comprehensive  incoming 
spectra  such  as  Pierson  -  Moskwitz  or  ionswap  (17).  More  realistic  «!>„  does  not  greatly  complicate  the 
procedure  The  basic  required  relations  have  been  derived  in  this  report 

£.  Short  Wavelength  Contribution  to  the  Surface  Energ) 

Shon  of  solving  for  the  complete  Fourier  spectral  distribution  throughout  the  interaction  zone  we 
can  also  cbrracterize  the  contributions  from  vanous  ranges  of  wavelengths  to  the  surface  energy.  This 
necessitates  establishing  modified  limits  of  integration  foi  the  energy  integral  previously  developed. 
That  is,  we  must  find  the  locus  of  constant  wavelength  within  the  interaction  zone  in  cither  the  initial 
or  physical  plane.  Equation  (47a)  or  (^7b)  may  then  be  used  to  calculate  tne  surface  energy  contribu¬ 
tion  from  the  wavelength  ranges  of  interest.  We  have  performed  this  calculation  in  the  initial  plan  uti¬ 
lizing  £q.  t47b).  In  general,  there  are  two  types  of  wavelength  contours  in  initial  space:  closed  and 
open.  TSe  major  differences  in  these  cases  arise  in  that  portion  of  frequency-direction  space  where  the 
waves  are  propagating  against  the  current.  The  evaluation  of  the  dimerisional  w-avelcngth  contour 
proceeds  from  Eqs.  (42)  and  (43).  These  can  be  combined  as 


'  k^i\  —  A  i)»)^ 


and  therefore 


Solving  for  sin  we  obtain 

if  I  , 

sm  0^  -  —rr - I - =  sm  6*.  t.^6) 

A  y„  <■>«•  I  Ct^ee 

This  equation  provides  the  value  of  9«.  which  produces  from  a  frequency  ct>e»  a  wave  number  of  A:  at  a 
position  y  where  the  local  cun-eni  is  given  by  A(y),  and  it  therefore  defines  the  required  locus  in 
the  initial  plane 

Equation  (56)  can  be  used  .o  establish  the  required  limits  of  integration  for  the  surface  energy  as 
discussed  above  Of  course,  the  limitations  imposed  by  Eq  (56)  hdva  to  be  incorporated  wah  respect 
to  those  limitations  already  associated  with  wave  breaking,  md  wave  reflection,  9^,  !n  thos,'  cases 
where  9^  from  Eq.  (56)  for  the  wave  number  boundary  are  intermediate  of  and  9^,  for  wave 
breaking  and  wave  reflection,  the  integration  is  from  9^  to  9^,  to  characienze  contnbutions  for  those 
wavelength  smaller  than  2ir/k  only.  This  sort  of  wave  length  decomposition  allows  a  partial  extraction 
of  the  spectral  densii)  distnbulion  The  full  procedures  established  above  for  the  numerical  evaluation 
of  the  spectral  density  will  provide  accurate  estimates  o“  the  Founer  coeRlcienis  of  the  free  surface 
elevation  These  modifications  are  now  in  progress  and  the  results  of  this  capability  will  be  reported 
soon.  The  remainder  of  this  report  will  concentrate  on  ihe  partial  spectral  knowledge  associated  with 
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the  k  contours  and  <I>„  (<o„,  fl„)  —  1.  Thai  is,  we  will  determine  the  surface  energy  contribuiions  for 
the  small  wavelength  range.  Subsequently,  these  computations  will  be  repeated  for  more  realistic 
incident  spectral  density  distributions. 

III.  RESULTS 

In  this  section  we  will  discuss  briefly  some  of  the  results  that  have  been  obtained  from  our  studies 
of  the  refraction  of  a  random  wave  field  by  a  simple  surface  current  shear  i/(y).  This  is  expected  to  be 
the  dominant  behavior  of  the  wake  velocity  field  of  a  surface  vessel.  These  results  have  bwn  selected 
to  provide  some  initial  estimates  for  the  effectiveness  of  this  wave-current  interaction  meciianism  for 
the  elimination  of  short  wavelength  components  in  the  downstream  region.  Short  wave  components 
may  have  their  origin  in  the  incoming  wave  spectrum  as  short  waves  or  they  may  be  general  from 
longer  vases  by  the  refraction  of  the  incoming  wave  field.  It  is  useful  to  distinguish  between  and  to 
quantify  these  two  origins. 

Figures  1  to  4  display  the  variation  of  mean-square  surface  elevation  energy  for  wavelengths  less 
than  Xr  (greater  than  as  a  function  of  position  within  different  strength  parabolic  wake  structures. 
iHOiXf  K,  -  1  -  k|  and  Xi  represent  wave  length  limits  for  the  uniform  spectra]  density  of 

the  incoming  waves  and  Xr.  represents  the  maximum  wavelength  (mininaum  wave  number 
—  2ir/Xr.)  which  is  mcluded  in  the  calculated  wave  field  energy,  £^.  presents  the  maxiroum  of 
the  current  velocity  profile.  It  is  ser  a  from  those  cases  where  Xj.  <"  Xj  <  Xi  or  Xj.  >  Xj>  Xj  that  the 
short  wavelength  components  may  be  crested  from  longer  wavelengths  Also  from  the  olner  cases  it  is 
seen  that  the  short  wavelength  components  originally  present  are  reduced  by  the  refraction  process. 

Figure  1  shows  a  comparison  for  fixed  values  of  Xc,  Xi  and  Xj  of  the  variation  of  the  root  mean 
square  surface  energy  for  wavelengths  less  than  Xj,  -  6.28  cm  produced  from  oncoming  wavelengths 
between  628  cm  and  12.56  cm  with  <I>«  -  i  for  wakes  with  various  maximum  velocities.  For  these 
cases,  the  short  wavelength  components  must  be  generated  by  the  refraction  process.  More  energetic 
short  vrave  components  are  generated  for  the  wake  with  a  maximum  velocity  of  *0  cm/sec  then  for.that 
with  a  maximum  of  30  cm/sec.  Both  the  30  and  40  cm/sec  maxima  are  not  strong  enough  to  produce 
any  wave  breaking  according  to  the  conservatSe  wave  breaking  criteria  given  by  Eqs.  (52a)  end  (52b). 
The  wakes  with  maximum  velocity  of  100  ap  .  200  cm/sec  do  cause  wave  breaking.  This  is  seen  for  the 
near  centerline  energy  which  is  reduced  in  those  cases  while  the  pioduction  of  short  waves  at  the 
current  structure  edges  is  higher.  Aiso  included  in  the  figure  is  the  variation  of  the  full  energy  over  all 
wavelengths  Ej-  foi  the  four  maximum  currents  cases  The  suongest  currents  produce  a  larger  reduc- 
uon  in  surface  mean  square  elevation  energy.  Note  that  in  those  cases  where  wave  breaking  occurs 
Ef! Else  1*®  than  unity  at  the  wake  boundaries,  ylyma  -  ±  1- 

The  same  energy  variation  with  position  in  the  current  is  shown  in  Fig.  2  for  the  cases  wheie  X/, 
IS  intermediate  between  Xi  and  Xj  and  for  two  values  of  maximum  current  both  of  which  do  not  induce 
wave  breaking.  The  balancing  effects  of  destruction  and  production  of  short  wave  length  components 
due  to  refraction  is  apparent  Figure  3  shows  these  effects  for  a  case  where  wave  breaking  is  present 
The  parameters  of  Fig.  3  are  Xi,  -  6.28  cm,  Xi  —  628  cm  and  Xi  -  6.28  cm.  The  maximum  current 
strength  is  —  70  cm/sec.  The  upper  curve  represents  the  variation  of  the  full  surface  energy  while 
the  lower  curve  represents  the  variation  of  the  surface  energy  for  wavelengths  less  than  6  2s  cm.  For 
the  majority  of  the  structure,  the  short  wave  ereigy  is  c-proximaiely  one-half  the  full  energy. 

In  Fig  4  we  display  the  same  variations  for  a  maximum  current  strength  of  ■■  20  cm.'sec  with 
X|  “  628  cm  and  Xj  -  Xr  -  6.28  cm.  No  breaking  occurs  for  this  .-ase.  The  short  wave  length  energy 
does  not  have  a  depression  near  the  maximum  velocity  at  the  centerline  of  the  current,  y  —  0. 


IV.  SUMMARY 


A  model  !>as  been  ouilined  in  this  report  for  characterizing  the  modification  of  an  incoming  ran¬ 
dom  ocean  wave  field  hy  a  surface  wake.  .At  the  most  fundamental  level,  we  can  compute,  for  a  white 
noise-like  incoming  spectrum  the  mean  square  surface  elevation  energy  and  its  variation  throughout 
the  wake  structure.  On  the  next  level  of  complication,  we  can  subAvide  the  .esults  into  various 
wavelength  regions  and  begin  to  describe  some  of  the  characteristics  of  die  spectral  changes.  Then, 
more  realistic  incoming  spectral  distributions  such  as  Pierson-Moskoiritz,  Jonswap,  etc.  can  be  incor¬ 
porated  into  the  model.  Lastly  and  most  importantly  we  will  be  able  to  determine  for  each  of  the  cases 
above  the  Founer  descitption  of  the  surface  elevation  field. 

From  the  representative  cases  which  have  been  considered  here,  it  is  suggested  that  a  significant 
reduction  in  the  ort  wavelength  components  of  the  incident  wave  field  is  possible  from  the  interac¬ 
tion  between  oce  j  surface  waves  and  currents  This  is  particularly  true  for  the  waves  with  components 
in  the  direction  of  the  current  Those  wave  components  oppositely  directed  to  the  current  are  reduced 
in  wavelength  by  the  interaction.  However,  they  are  refracted  so  that  they  are  turned  away  from  the 
current  direction  This  win  show  mere  clearly  when  we  are  able  'o  determine  the  relative  amplitudes  of 
the  components  of  the  frequency-direction  spectral  density.  The  relative  absence  of  components  with 

angles  close  to  0  — -y  (opposite  to  current )  may  be  important  in  analyzing  the  response  of  a  particu¬ 

lar  remote  sensing  sj’stem.  This  directional  distribution  may  be  more  important  than  the  relative  popu¬ 
lation  of  short  wave  lengths. 
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Fig.  1  —  Variation  of  nonnalized  mean  square  surface  displacement  energy  as  a  function  of  position  in 
a  parabolic  wake  structure  for  different  strength  wakes,  y„.  represents  the  surface  energy  for 
wavelengths  less  than  Xi-  £r  represenls  the  surface  energy  for  all  wavelengths  considered  and  Ejk 
represents  the  total  surface  energy  for  the  incoming  spectrum.  Incoming  spectrum  -  4>»  -  1  tor 
12.56  cm  <  X  <  628  cm. 


